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Abstract 

The ability to perform high-precision one- and two-qubit operations is 
sufficient for universal quantum computation. For the Loss-DiVincenzo 
proposal to use single electron spins confined to quantum dots as qubits, 
it is therefore sufficient to analyze only single- and coupled double-dot 
structures, since the strong Heisenberg exchange coupling between spins 
in this proposal falls off exponentially with distance and long-ranged dipo- 
lar coupling mechanisms can be made significantly weaker. This scal- 
ability of the Loss-DiVincenzo design is both a practical necessity for 
eventual applications of multi-qubit quantum computing and a great con- 
ceptual advantage, making analysis of the relevant components relatively 
transparent and systematic. We review the Loss-DiVincenzo proposal for 
quantum-dot-confined electron spin qubits, and survey the current state of 
experiment and theory regarding the relevant single- and double- quantum 
dots, with a brief look at some related alternative schemes for quantum 
computing. 
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1 INTRODUCTION 

Recent advances in semiconductor spintronics and, more specifically, spin-based 
quantum computing in solid-state systems, have encouraged significant research 
efforts in the last years (Prinz 1998, Wolf et al. 2001, Awschalom et al. 2002). 
Much of this research is motivated by pressure on the electronics industry 
to maintain Moore's-law growth in systems with components that are very 
quickly approaching the nanoscale, where quantum mechanics becomes impor- 
tant (ITRS 2005). Additionally, nanoscale devices provide a unique opportunity 
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to study the fundamental physics of quantum phenomena in a controllable en- 
vironment. 

Independent of the particular motivation, if quantum information processing 
is to progress beyond basic proof-of-principle experiments, it must be based on 
a viable, scalable qubit (a quantum mechanical two-level system, which can be 
placed in an arbitrary superposition of its basis states: \ip) = a |0) + b |1)). The 
two states of single electron spins (||) = |0) and ||) = |1)), confined to semicon- 
ductor quantum dots (the Loss-DiVincenzo proposal), are one such qubit (Loss 
and DiVincenzo 1998). These qubits are viable, in the sense that they make use 
of fabrication techniques and electrical control concepts that have been devel- 
oped over the last five decades in research laboratories and industry. The secret 
to scalability in the Loss-DiVincenzo proposal lies in local gating; this proposal 
would implement gating operations through the exchange interaction, which 
can be tuned locally with exponential precision, allowing pairs of neighboring 
qubits to be coupled and decoupled independently This is to be contrasted with 
proposals that make use of long-ranged interactions (e.g., dipolar coupling) for 
which scalability may be called into question. The local, tunable nature of inter- 
qubit interactions in the Loss-DiVincenzo proposal is what makes it possible to 
consider first isolated one-qubit (single-quantum-dot), then isolated two-qubit 
(doublc-quantum-dot) systems. Once single- and double- quantum dots are 
understood, along with environmental couping mechanisms, a quantum compu- 
tation can proceed through a scries of one- and two-qubit operations, without 
great concern regarding interactions between three, four, and more qubits. 

There are many other proposals for qubits and associated quantum control 
processes. Some examples include various proposals that use superconducting 
devices (for reviews, see (Makhlin et al. 2001, Burkard 2004)), proposals for "adi- 
abatic quantum computing", in which quantum computations are performed 
through adiabatic manipulation of coupling constants in physically realizable 
Hamiltonians (Farhi et al. 2000, Farhi et al. 2001, Wu et al. 2005) (which might 
be used to perform fast quantum simulations of, e.g., superconducting pair- 
ing models (Wu et al. 2002)), electron spin qubits encoded in two-spin states 
(Levy 2002) or many-spin chains (Meier et al. 20036, Meier et al. 2003a) (recent 
work showing that such spin chains can be built-up atom-by-atom on a surface 
(Hirjibehedin et al. 2006) is a possible first step to implementing such a pro- 
posal), cavity-QED schemes (Sleator and Weinfurter 1995, Domokos et al. 1995), 
trapped- ion proposals (Cirac and Zoller 1995), etc. Each of these proposals has 
advantages and disadvantages. Here we do not compare the relative merits of 
all proposals, but instead focus on proposals involving electron spins confined 
to quantum dots. 

Before a quantum computation can begin, the qubits in a working quantum 
computer must be initialized to some state, e.g. |0). These qubits must be 
sufficiently isolated from the surrounding environment to reduce decoherence, 
there must be some way to perform fast single- and two-qubit operations in 
a time scale much less than the qubit decoherence time, and it must be pos- 
sible to read out the final state of the qubits after any quantum computation 
(DiVincenzo 2000). In the following sections, we address these issues and others 
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which arc important for the Loss-DiVincenzo proposal. In the process, we sur- 
vey some recent work on quantum computing with electron spins in quantum 
dots. 

Due to the rapid development of this field, there have been many recent 
reviews on quantum-dot quantum computing. Among these numerous reviews, 
there has been work that focuses on single-electron charge qubits in double 
dots (Fujisawa et al. 2006), the implementation of single-electron spin reso- 
nance (ESR), and the molecular wavefunctions of coupled double dots (van der 
Wiel ct al. 2006), various proposals for spin-based quantum computing (Cerletti 
et al. 2005), silicon-based proposals for quantum computing (Koiller et al. 2005), 
general quantum computing in the solid state, including both quantum dots 
and superconducting systems (Burkard 2004), experiments and experimental 
proposals for quantum-dot-confined electrons (Engel et al. 2004 a), the many 
coupling schemes and decoherence mechanisms for quantum-dot spin qubits 
(Hu 2004), and optical properties of quantum dots (Hohencster 2004). In this 
review, we analyze the Loss-DiVincenzo proposal from the viewpoint that this 
proposal can be decomposed into first single and then double quantum dots, 
with a special emphasis on double-dot physics. 

This review is organized as follows: in Section [3 we give a brief summary of 
the Loss-DiVincenzo proposal for quantum computing. In Section [21 we discuss 
the characterization and manipulation of spin and charge states of electrons in 
single quantum dots. Section 0] contains a description of double quantum dots 
that emphasizes the single-electron regime, which is relevant for quantum-dot 
quantum computing. In Scction0we survey important decoherence mechanisms 
for electron spins in single- and double- quantum dots. In Sectional we briefly 
review some proposals for the generation and detection of nonlocal entanglement 
of electron spins in nanostructures, and in Section [7| we conclude with a brief 
summary of important topics for future study. 



2 SPINS IN QUANTUM DOTS: AN OVERVIEW 
OF THE LOSS-DIVINCENZO PROPOSAL 

In the original Loss-DiVincenzo proposal, the qubits are stored in the two spin 
states of single confined electrons. The considerations discussed in (Loss and 
DiVincenzo 1998) are generally applicable to electrons confined to any structure 
(e.g. atoms, molecules, defects, etc.), although the original proposal focused 
on applications in gated semiconductor quantum dots, as shown in Figure ^ 
Voltages applied to the top gates of such structures provide a confining potential 
for electrons in a two-dimensional electron gas (2DEG), below the surface. A 
negative voltage applied to a back-gate depletes the 2DEG locally, allowing 
the number of electrons in each dot to be reduced down to one (the single- 
electron regime). Advances in materials fabrication and gating techniques have 
now allowed for the realization of single electrons in single vertical (Tarucha 
et al. 1996) and gated lateral (Ciorga et al. 2000) dots, as well as double dots 
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Figure 1: A double quantum dot. Top-gates are set to an electrostatic voltage 
configuration that confines electrons in the two-dimensional electron gas (2DEG) 
below to the circular regions shown. Applying a negative voltage to the back- 
gate, the dots can be depleted until they each contain only one single electron, 
each with an associated spin-1/2 operator Sl(r) for the electron in the left 
(right) dot. The |f) and |j) spin-1/2 states of each electron provide a qubit 
(two-level quantum system). 

(Elzerman et al. 2003, Hayashi et al. 2003, Petta et al. 2004). 

Initialization of all qubits in the quantum computer to the Zeeman ground 
state |t) = |0) could be achieved by allowing all spins to reach thermal equi- 
librium at temperature T in the presence of a strong magnetic field B, such 
that |<7/is-B| > ksT, with ^-factor g < 0, Bohr magneton /xb, and Boltzmann's 
constant fcs (Loss and DiVincenzo 1998). For further initialization schemes, see 
Section 13.21 below. 

Once the qubits have been initialized to some state, they should remain in 
that state until a computation can be executed. In the absence of environmental 
coupling, the spins-1/2 of single electrons are intrinsic two-level systems, which 
cannot "leak" into higher excited states. Additionally, since electron spins only 
couple to charge degrees of freedom indirectly through the spin-orbit (or hyper- 
fine) interactions, they are relatively immune to fluctuations in the surrounding 
electronic environment. 

Single-qubit operations in the Loss-DiVincenzo quantum computer could be 
carried out by varying the Zeeman splitting on each dot individually (Loss and 
DiVincenzo 1998). It may be possible to do this through g-factor modulation 
(Salis et al. 2001), the inclusion of magnetic layers (Myers et al. 2005) (see also 
Figure 0), modification of the local Overhauser field due to hypcrfinc couplings 
(Burkard et al. 1999), or with nearby ferromagnetic dots (Loss and DiVincenzo 
1998). There are a number of alternate methods that could be used to perform 
single-qubit rotations (see Section l3"2)l . 

Two-qubit operations would be performed within the Loss-DiVincenzo pro- 
posal by pulsing the exchange coupling between two neighboring qubit spins 
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Figure 2: A series of exchange-coupled electron spins. Single-qubit operations 
could be performed in such a structure using electron spin resonance (ESR), 
which would require an rf transverse magnetic field B^, and a site-selective 
Zeeman splitting g(x)fiBB±, which might be achieved through ^-factor modu- 
lation or magnetic layers. Two-qubit operations would be performed by bringing 
two electrons into contact, introducing a nonzero wavefunction overlap and cor- 
responding exchange coupling for some time (two electrons on the right). In 
the idle state, the electrons can be separated, eliminating the overlap and cor- 
responding exchange coupling with exponential accuracy (two electrons on the 



"on" to a non-zero value (J(t) = Jo 7^ 0, t G {— r s /2 . . . t s /2}) for a switching 
time t s , then switching it "off" (J(t) = 0, t ^ {— r s /2 . . .r s /2}). This switching 
can be achieved by briefly lowering a center-gate barrier between neighboring 
electrons, resulting in an appreciable overlap of the electron wavefunctions (Loss 
and DiVincenzo 1998), or alternatively, by pulsing the relative back-gate voltage 
of neighboring dots (Petta et al. 2005 a) (see Section PPf . Under such an opera- 
tion (and in the absence of Zeeman or weaker spin-orbit or dipolar interactions), 
the effective two-spin Hamiltonian takes the form of an isotropic Heisenberg ex- 
change term, given by (Loss and DiVincenzo 1998, Burkard et al. 1999) 



where Sl(r) is the spin-1/2 operator for the electron in the left (right) dot, as 
shown in Figure ^ The Hamiltonian H ex (t) generates the unitary evolution 
U(4>) — exp [— «0Sl • Sr], where = J J(t)dt/h. If the exchange is switched 
such that (f> = J J(t)dt/H = Jot s /H = tt, U((f>) exchanges the states of the two 
neighboring spins, i.e.: U(ir) |n,n') = |n',n), where n and n' are two arbitrar- 
ily oriented unit vectors and |n, n') indicates a simultaneous eigenstate of the 
two operators Sl ■ n and Sr. • n'. U(ir) implements the so-called SWAP opera- 
tion. If the exchange is pulsed on for the shorter time r s /2, the resulting op- 
oration U(tt/2) = (U(n)) ' is known as the "square-root-of-SWAP" (vswap). 
The v / SWAP operation in combination with arbitrary single-qubit operations 



left). 



H ex (t) - J(t)S h ■ S R , 



(1) 



5 



3 SINGLE QUANTUM DOTS 



is sufficient for universal quantum computation (Barenco et al. 1995a, Loss 
and DiVincenzo 1998). The \/swap operation has now been successfully im- 
plemented in experiments involving two electrons confined to two neighboring 
quantum dots (as in Figure 0) (Petta et al. 2005a, Laird et al. 2005). Errors 
during the v / SWAP operation have been investigated due to nonadiabatic tran- 
sitions to higher orbital states (Schliemann et al. 2001, Requist et al. 2005), 
spin-orbit-interaction (Bonesteel et al. 2001, Burkard and Loss 2002, Stepa- 
nenko et al. 2003), and hyperfine coupling to surrounding nuclear spins (Petta 
et al. 2005a, Coish and Loss 2005, Klauser et al. 2005, Taylor et al. 2006). The 
isotropic form of the exchange interaction given in Equation Q is not always 
valid. In realistic systems, a finite spin-orbit interaction leads to anisotropic 
terms which may cause additional errors, but could also be used to perform 
universal quantum computing with two-spin encoded qubits, in the absence of 
single-spin rotations (Bonesteel et al. 2001, Lidar and Wu 2002, Stepanenko and 
Bonesteel 2004, Chutia et al. 2006) (see also Section El below). 

In the Loss-DiVincenzo proposal, readout could be performed using spin- 
to-charge conversion. This could be accomplished with a "spin filter" (spin- 
selective tunneling) to leads or a neighboring dot, coupled with single-electron 
charge detection (see also Section l3~^l below). 

3 SINGLE QUANTUM DOTS 

The fundamental element of information in a quantum computer is the quantum 
bit, or qubit. The qubits of the Loss-DiVincenzo proposal (Loss and DiVincenzo 
1998) are encoded in the two spin states ||) and ||) of single electrons confined 
to quantum dots. 

There are many different types of quantum dot that can be manufactured, 
each with distinct characteristics. Gated lateral quantum dots (as shown in 
Figures ^ an d El offer the benefit that their shape and size can be controlled 
to suit a particular study, and the tunnel coupling between pairs of these dots 
can be tuned in a straightforward manner: by raising or lowering the barrier 
between the dots. Gated vertical dots (Tarucha et al. 1996) are created by 
etching surrounding material to form a pillar structure, with vertical confine- 
ment provided by a double-barrier heterostructure. Vertical dots allow for the 
controlled fabrication of quantum dots with large level spacing, although tun- 
ability of the coupling in these structures is restricted due to the fabrication 
process. To resolve this issue, hybrid laterally-coupled vertical double quantum 
dots have been manufactured, in which the inter-dot tunnel coupling is control- 
lable (Hatano et al. 2005). Self-assembled quantum dots are yet another type of 
dot that can be used for quantum information processing. Self-assembled dots 
form spontaneously during epitaxial growth due to a lattice mismatch between 
the dot and substrate materials. These dots can be made with very large single- 
particle level spacings, but typically form at random locations, which makes 
controlled coupling through a tunnel junction difficult. Such dots can, how- 
ever, potentially be coupled with optical cavity modes (Imamoglu et al. 1999), 
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and new techniques have now allowed the fabrication of cavities with modes that 
couple maximally directly at the positions of isolated dots (Badolato et al. 2005). 

In the rest of this section we focus on lateral quantum dots, as shown in 
Figures^and|5| After a brief review of charge and spin control in single quantum 
dots, we will address issues specific to double quantum dots in Section 0] 

3.1 Charge control: Coulomb blockade 

To ensure a single two- level system is available to be used as a qubit, it is practi- 
cal to consider single isolated electron spins (with intrinsic spin 1/2) confined to 
single orbital levels. A natural first step to implementing the Loss-DiVincenzo 
proposal was therefore to demonstrate control over the charging of a quantum 
dot electron- by-electron in a single gated quantum dot. This is typically done 
by operating a quantum dot in the Coulomb-blockade regime, where the energy 
for the addition of an electron to the quantum dot is larger than the energy 
that can be supplied by electrons in the source or drain leads. In this case, the 
charge on the quantum dot is conserved, and no electrons can tunnel onto or off 
of the dot. For a general review of Coulomb blockade phenomena and the char- 
acterization of many-electron states in single quantum dots, see (Kouwenhoven 
et al. 2001). 

3.2 Spin control: Initialization, operations, and readout. 

As mentioned in Section [21 initialization of all electron spins to the "up" state 
|t) could be achieved by allowing all spins to equilibrate in a strong magnetic 
field. Depending on the particular architechture, this may take a long time 
or it may be inconvenient to have large magnetic fields in the region of the 
apparatus. Initialization could also be achieved through spin-injection from a 
ferromagnet, as has been performed in bulk semiconductors (Fiederling et al. 
1999, Ohno et al. 1999), with a spin-polarized current from a spin-filter device 
(Prinz and Hathaway 1995, Prinz 1998, Loss and DiVincenzo 1998, DiVincenzo 
1999, Recher et al. 2000), or by optical pumping (Cortez et al. 2002, Shabaev 
et al. 2003, Gywat et al. 2004, Bracker et al. 2005), which has now allowed the 
preparation of spin states with very high fidelity, in one case as high as 99.8% 
(Atature et al. 2006). 

Single-qubit operations can be performed in the Loss-DiVincenzo proposal 
whenever the Zeeman energy of the quantum-dot spins can be tuned locally, 
as mentioned in Section Alternative single-qubit-rotation schemes may re- 
quire global magnetic field gradients (Wu et al. 2004, Tokura et al. 2006), ESR 
(see Figure |2J) or, in the presence of spin-orbit interaction, electric-dipole spin 
resonance (EDSR) techniques. EDSR has been analyzed in great detail for 
two-dimensional systems in theory (Rashba and Efros 2003, Duckheim and 
Loss 2006) and experiment (Kato et al. 2004), and can also be applied to 
lower-dimensional systems (quantum wires and quantum dots) (Levitov and 
Rashba 2003, Golovach et al. 2006), with the advantage that single-qubit opera- 
tions could then be performed using fast all-electrical control. New experiments 
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have now shown that it may be possible in practice to perform single-spin op- 
erations on single quantum dots using ESR, as depicted in Figure [21 (Koppens 
et al. 2006). 

As mentioned in Section |2 quantum-dot spin readout can be performed 
using a spin filter. Experimentally, spin filters have been reported in the open 
(Potok et al. 2002) and Coulomb-blockade regimes (Folk et al. 2003), and have 
even been used to determine the longitudinal spin decay (T±) time (Hanson 
et al. 2003, Hanson et al. 2004) using an n-shot readout scheme, which has been 
analyzed in detail (Engel et al. 2004 6). A single- shot readout has also been 
demonstrated (Elzerman et al. 2004) and improved upon (Hanson et al. 2005). 
Non-invasive readout schemes using spin-to-charge conversion and quantum- 
point-contact (QPC) measurements have been used on two-spin encoded qubits 
(Johnson et al. 2005a, Petta et al. 2005 a, Petta et al. 20056, Johnson et al. 
20056). 

To measure the transverse spin coherence time T2, there have been proposals 
to perform ESR and detect the resulting resonance in stationary current (Engel 
and Loss 2001), changes in the resistivity of a neighboring field-effect transistor 
(FET) (Martin et al. 2003), optically (Gywat et al. 2004), or from current noise 
(Schaefers and Strunz 2005). ESR in single quantum dots has not yet been 
observed, in part because it is challenging to generate high-frequency magnetic 
fields with sufficient power for single-spin manipulation without "heating" elec- 
trons on the quantum dot or in the surrounding leads through the associated 
electric field (van der Wiel et al. 2006). Recent experiments that employ a dou- 
ble quantum dot in the spin-blockade regime may have overcome this problem 
(Koppens et al. 2006) (see also the discussion on spin blockade near the end of 
Section El below). 

4 DOUBLE QUANTUM DOTS 

Single qubits are the fundamental unit of quantum information in quantum 
computing. However, universal quantum computation still requires both single- 
qubit and two-qubit operations (Barenco et al. 19956). In the Loss-DiVincenzo 
proposal, two-qubit gates are performed with exchange-coupled electron spins 
confined to two neighboring quantum dots (double dots). Double dots are also 
important for encoded qubits (Levy 2002), in which qubits are encoded into 
a two-dimensional pseudospin-1/2 subspace of a four-dimensional two-electron 
spin system. 

In this section we discuss characterization and manipulation techniques that 
are commonly used to extract microscopic parameters of double quantum dots. 
In Section 14.11 we review the charge stability diagram, and illustrate its con- 
nection to a commonly used microscopic model Hamiltonian. In Section 14.21 
we review work on the coherent coupling of double quantum dots, which is 
required to generate a large exchange interaction for two-qubit gating. In Sec- 
tion 14.31 we discuss the use of double quantum dots as two-qubit gates, and 
in Section 14.41 we review some work on using double quantum dots to control 
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single "encoded" qubits (Levy 2002), a topic which has now come into vogue 
(Petta et al. 2005a, Taylor et al. 2005, Burkard and Imamoglu 2006, Hanson 
and Burkard 2006). 

4.1 The double-dot charge stability diagram 

Just as transport through a single quantum dot and Coulomb blockade phe- 
nomena give information about the orbital level spacing, charging energy, and 
spin states of single quantum dots, similar studies can be carried-out on dou- 
ble quantum dots. Whereas for single dots, transport phenomena are typically 
understood in terms of one-dimensional plots of conductance versus gate volt- 
age, the primary tool used to understand double quantum dots is the double-dot 
charge stability diagram. The stability diagram is a two-dimensional plot of cur- 
rent or differential conductance through the double dot or through a neighboring 
QPC, given as a function of two independent back-gate voltages (one applied 
locally to each dot) . The plot differentiates regions where the double-dot ground 
state has a charge configuration (N\,N2), for various Ni,N2, where Ni is the 
number of charges on the left dot and N2 is the number of charges on the right. 
Transport through double quantum dots and the relevant charge stability dia- 
gram has been discussed thoroughly in (van der Wiel et al. 2003). In the rest of 
this section, we review some features of the double-dot stability diagram with 
an emphasis on the connection to a model Hamiltonian that is commonly used 
in the literature (Klimeck et al. 1994, Pals and MacKinnon 1996, Golden and 
Halperin 1996, Ziegler et al. 2000). 

An isolated double quantum dot is described by the Hamiltonian 

H dd = Hc + Ht + Hs, (2) 

where Hq gives the single-particle and inter-particle charging energies as well 
as the orbital energy, Ht is the inter-dot tunneling term due to a finite over- 
lap of dot-localized single-particle wavefunctions, which ultimately gives rise to 
exchange, and iJg contains explicitly spin-dependent terms, which may include 
spin-orbit interaction, dipolc-dipole interaction, and the contact hypcrfine inter- 
action between the confined electron spins and nuclear spins in the surrounding 
lattice. 

There are several approaches that can be taken to writing the various com- 
ponents of the double-dot Hamiltonian H^d, corresponding to several degrees 
of microscopic detail. In the simplest form, the Hubbard model, details of the 
electron wavefunctions are neglected and the Coulomb interaction is given only 
in terms of on-site and nearest-neighbor terms. Since this description relies 
only on very few parameters, it is the most commonly used in the literature 
on transport phenomena through quantum dots. The shape of the confining 
potential, quantum-dot localized wavefunctions, and form of the Coulomb in- 
teraction may become important in certain circumstances, in which case it is 
more appropriate to apply either the Heitler-London method (which neglects 
doubly-occupied dot levels), or the Hund-Mulliken method, which includes the 
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Figure 3: Ground-state configuration for a double quantum dot with large or- 
bital and charging energies, and negligible dot-lead and interdot coupling. /i s (d) 
is the source (drain) chemical potential, Vin) is the left (right) local dot poten- 
tial, which is related to applied gate potentials by a linear transformation (see 
Equation (JSJ , below) , and both dots are assumed to have the same uniform level 
spacing Kujq. 

effects of double-occupancy. These methods predict, for instance, a variation 
of the interdot exchange interaction through zero with increasing out-of-planc 
magnetic field (Burkard et al. 1999). Experimentally, it has been confirmed 
that the exchange coupling can be tuned with an out-of-plane magnetic field in 
single vertical (Fujisawa et al. 2002) and single lateral quantum dots (Zumbiihl 
et al. 2004), which behave effectively as double-dot structures. Here we ignore 
these effects and focus on the simplest Hubbard model that reproduces much of 
the double-dot physics that can be seen in transport phenomena. 

We model the Coulomb interaction with simple on-site (J7i(2) for the left 
(right) dot) and nearest-neighbor (U 1 ) repulsion. The single-particle charging 
energy is given in terms of a local dot potential Vin) ■ The charging Hamiltonian 
is then 

hc = \y< UiNi w - + U ' N ^ - m ViNi + eiknik > ( 3 ) 

l l kl 

where Ni = nik counts the total number of electrons in dot Z, with nik = 
^2 a dj k(T dika, and here dik a annihilates an electron on dot I, in orbital fc, with 
spin a. eik is the energy of single-particle orbital level fc in dot I, which gives 
rise to the typical orbital level spacing eik+i — (ik ~ ^>o (see Figure[3j). 

Within the capacitive charging model described by the equivalent circuit 
in the inset of Figure EJa), the microscopic charging energies are related to 
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Figure 4: Stability diagram plotted in terms of (a) local dot potentials V\ t 2 and 
(b) applied gate potentials V g \^, with on-site charging energies Ui = U, I = 1, 2, 
nearest-neighbor charging energy U', and dot orbital level spacing fuvo satisfying 
U : hujQ : U' = 3 : 2 : 1. In addition, for (b) we have assumed the voltage scaling 
factors are the same for both dots, and are given by a\ = ct2 = a = 1/2. 
(a) inset: capacitive charging model for a double quantum dot, indicating the 
source (drain) chemical potential /i s (d), the charge on the left (right) dot Qu 2)1 
the capacitances to source (drain) C s (<j), the mutual capacitance C m , and gate 
capacitances C g i, 2 . (b) Horizontal lines in the |e| V^( 2 ) plane become skewed 
with slope SV g i/SV g2 = — CiCg 2 /C m Cgi when plotted versus |e| ^1(2)- 
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capacitances by (Ziegler et al. 2000, van der Wiel et al. 2003) 

C1C2 e 2 2e 2 C m 

where Ci = C s + C m + C g i , C2 = Cd + C m + C g 2 , and all capacitances are denned 
in the inset of Figure^a). In experiments, the local quantum dot potentials Vi t 2 
are controlled indirectly in terms of gate voltages Vgi.2, which are capacitively 
coupled to the dots through gate capacitances C g i,2- For fixed quantum-dot 
charges (Qi,Q2) = — \e\ (N\,N2) = const., differences in the dot voltages AVi 
and AV2 are related to differences in the gate voltages AVgi and AVg2 through 
(Ziegler et al. 2000, van der Wiel et al. 2003) 



Ci -C m \(AV 1 \ = ( C gl AV gl 
-C m C 2 ) \ AV 2 J \ C e2 AV g2 



(5) 



The double-dot stability diagram can then be given equivalently as a two- 
dimensional plot with energy axes |e|Vi, |e | V2 , or with axes |e|Vgi, |e|Vg2, 
which are skewed and stretched with respect to the original axes according to 
the transformation given in Equation The end effect is that parallel horizon- 
tal (vertical) lines in the |e| Vi(2) plane separated by a distance <LV2(x) transform 
to skewed parallel lines, separated by dV g 2(i) = dV 2 (i) / ce 2 (i) along the horizontal 
(vertical) of the new coordinate system, where (see Figure : 

<*i = % 1 = 1,2. (6) 

Additionally, horizontal lines in the |e| ^(2) plane become skewed with a slope 
<5Vgi/(5Vg2 = — C m Cg2/C2Cgi (see Figure E^b)), and vertical lines are skewed 
with slope (5Vgi/(5Vg2 = — CiC g 2/C m C g i. 

The Hamiltonian in Equation @ conserves the number of electrons on each 
dot: [Hc,Ni] — 0, so we label the ground state by the two dot occupation 
numbers, (JVijJVg), and indicate where each configuration is the ground state 
in Figure 21 for equivalent quantum dots that satisfy ol\ = 012 = a = 1/2, 
U\ = U2 = U, eik+i — £;fc = fo^o f° r a U and U : Hljq : U' = 3 : 2 : 
1. The charge stability diagram produces a "honeycomb" of hexagons with 
dimensions that are determined by three typical energy scales: (1) The on- 
site replusion U, (2) the nearest-neighbor repulsion U' , and (3) the typical 
orbital energy Hujq. Figure 0| assumes a ground-state electron filling as shown 
in Figure |3 with constant orbital energy Hujq. In this case, the orbital energy 
appears in the dimensions of only every second honeycomb cell of the stability 
diagram, along the horizontal or vertical direction, since the spin-degenerate 
orbital states fill with two electrons at a time according to the Pauli principle. 
This even-odd behavior may not be visible in dots of high symmetry, where the 
orbital levels are manifold degenerate. Alternatively, the absence of an even-odd 
effect in low-symmetry single dots has previously been attributed to the absence 
of spin degeneracy due to many-body effects (Stewart et al. 1997, Fujisawa 
et al. 2001, van der Wiel et al. 2003). 
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Each vertex of a honeycomb cell corresponds to a triple-point, where three 
double-dot charge states are simultaneously degenerate. For a double dot con- 
nected to source and drain leads at low temperature, and in the absence of 
relaxation or photo-assisted tunneling processes, it is only at these points where 
resonant sequential transport can occur, through shuttling processes of the form 
(0,0) — ■> (1,0) — » (0,1) — > (0,0). This picture changes when a strong inter-dot 
tunnel coupling Ht is considered in addition. 

4.2 Molecular states in double dots 

Molecule-like states have been observed and studied in detail in two-electron 
single vertical (Fujisawa et al. 2002) and lateral quantum dots (Zumbiihl et al. 
2004) (the latter behave as an effective double-dot structure, showing good 
agreement with theory (Golovach and Loss 2004)). Evidence of molecular states 
forming in double quantum dots due to a strong inter-dot tunnel-coupling has 
also been found in a variety of systems (Schmidt et al. 1997, Schedelbeck et al. 
1997, Blick et al. 1998, Brodsky et al. 2000, Bayer et al. 2001, Ota et al. 2005, 
Hiittel et al. 2005, Fasth et al. 2005, Mason et al. 2004, Biercuk et al. 2005, 
Graeber et al. 2006). For example, molecular states have been observed in 
many-electron gated quantum dots in linear transport (Blick et al. 1998) (solid 
lines of Figure Elb)) and transport through excited states (Hiittel et al. 2005) 
(dashed lines in Figure [3b)). I n addition, molecular states have been observed 
in vertical-lateral gated double quantum dots (Hatano et al. 2005), gated dots 
formed in quantum wires (Fasth et al. 2005) and gated carbon-nanotube double 
dots (Mason et al. 2004, Biercuk et al. 2005, Graeber et al. 2006). A large inter- 
dot tunnel coupling is essential for generating a large exchange interaction J, 
and is therefore very important for the implementation of fast two-qubit gates 
in the Loss-DiVincenzo proposal. 

In this section, we analyze changes to the double-dot stability diagram that 
occur due to the inter-dot tunneling term H^. We focus on the relevant regime 
for quantum computing, where only a single orbital state is available for occu- 
pation on each quantum dot (the lower- left region of Figures 0Ja,b)). In the 
subspace of these lowest dot orbital states, Ht is given by: 



where t\2 is the tunneling amplitude between the two dots, and di a , I = 1,2, 
annihilates an electron in the lowest single-particle orbital state localized on 
quantum dot I with spin a. 

When the double dot is occupied by only N = 0, 1 electrons and is coupled 
weakly to leads, an explicit expression can be found for the current passing 
through a sequentially-coupled double dot, as shown in Figure OJa) (Ziegler 
et al. 2000, Graeber et al. 2006). It is straightforward to diagonalize He + Ht 
in the subspace of N — 1 electrons on the quantum dot. This gives the (spin- 




(7) 



a 
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Figure 5: (a) A tunnel-coupled double quantum dot, with tunneling amplitude 
ti2- The source and drain leads, at chemical potentials fi s and /id, are connected 
to the left and right dots through tunnel barriers with tunneling amplitudes t s 
and td, respectively. The left and right dots are set to local potentials V\ and 
V2 ■ (b) Modification of the stability diagram in the case of a significant tunnel 
coupling ti2 . To generate this figure we have chosen the ratio of tunnel coupling 
to the mutual (nearest-neighbor) charging energy to be t\ 2 jU' » 1/5. At solid 
lines, transport occurs via the double-dot ground state \E + ) and at dashed 
lines, additional transport can occur through the first excited state \E-) (sec 
Equations © and ffT7Tf> below). 

degenerate) eigenenergies and corresponding eigenvectors: 

E±(A,e) = -±(A±Je* + 2t* 2 ), (8) 
\E±) = cosf%) |l,0)+sin(^) |0,1> , (9) 




Ht) = V^^Hv^J- (10) 

Here, E±(A,e) is written in terms of new energy coordinates e, A, which are 
related to the old (voltage) coordinates through a rotation of the axes by 45° 
(see also Figure Efb)): 



A \ J_f 1 1W lelVr 

e J V2 V - 1 1 J \ 



(ii) 



We then define double-dot chemical potentials: 

H±(A,e)=E±(A,e)-E , (12) 

where Eq — is the energy of the (0, 0) charge configuration. In the pres- 
ence of a strong tunnel coupling, the eigenstates of the double dot are no 
longer labeled separately by the quantum numbers Ni, Instead, the sum 
N = Ni + N 2 is conserved. If we add to -ffdd the double-dot-lead coupling 
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Hamiltonian i?dd-£ = Y^k(7 taC ska^ + *dC dfc<r d2<7 + H-c, where c^ (d)fco . creates 
an electron in the source (drain), in orbital k with spin <r, then N can fluctu- 
ate between 1 and if the double-dot and lead chemical potentials are equal. 
We identify double-dot sequential tunneling processes as those that change the 
total charge on the double dot by one: N — > N ± 1 (Golovach and Loss 2004). 
One can evaluate golden-rule rates for all sequential-tunneling processes, taking 
the dot-lead coupling Hm-l as a perturbation to obtain the stationary current 
from a standard Pauli master equation (the Pauli master equation is valid for 
sufficiently high temperature, ksT > r s / d ), so that off-diagonal elements can be 
ignored in the double-dot density matrix). For weak dot-lead coupling, at low 
temperature k^T < Hljo, and at zero bias (fx = fi s = /id + A/i, with A/i — > 0), 
transport occurs only through the N — 1 ground state, with chemical potential 
/Lt_|_. The differential conductance near the N = 0, 1 boundary is then given by 

£ sin 2 (fl + )r s r d 

where f(E) = 1/ 1 + cxp ( J is the Fermi function at chemical potential /z 

and temperature T, f'(E) = df(E)/dE, and r s ( d ) = 2E£ |t s ( d )| 2 is the tunneling 
rate to the source (drain) with a lead density of states per spin v at the Fermi 
energy. If spin degeneracy is lifted, the quantity in brackets in Equation Ijl3|l is 
replaced by the familiar term — /'(/i + ) = 1/ 4fceTcosh 2 f 2fesT ) (Beenakker 
1991). The differential conductance (Equation 1)131) ) reaches a maximum near 
the point where the double-dot chemical potential matches the lead chemical 
potential, /i+(A, e) = which we indicate with a solid line in Figure |5jb). 
Transport through the excited state can occur where /x_(A,e) = /i, and when 
the bias A/i = \x s — fid or temperature T are sufficiently large to generate a 
significant population in the excited state |-EL). Dashed lines indicate where 
/U_(A, e) = /i in Figure 03Tb)- 

There are several qualitative changes to the double-dot stability diagram 
that take place in the presence of strong tunnel coupling. First, the number 
of electrons on each dot is not conserved individually. Instead, the sum N — 
Ni + N2 is conserved, which means that there are no longer lines separating, for 
example, the (1,0) and (0,1) states in Figure|Sfb). Second, sequential-tunneling 
processes allow current to be transported through the double-dot along the 
length of the "wings" that define the boundaries between TV and N ± 1-electron 
ground states. This is in contrast to the case where tn is weak, in which 
resonant sequential transport can only occur at triple points, where the shuttling 
processes of the type (0,0) — > (1,0) — > (0,1) — > (0,0) are allowed by energy 
conservation. 



4.3 Double dots as two-qubit gates 

The v / SWAP operation described in Section [21 requires significant control of the 
exchange coupling J. The value of J can be controlled by raising/lowering the 
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Figure 6: Energy-level spectrum for two electrons in a double quantum dot. 



inter-dot barrier, thus changing the tunnel coupling t\2 (Loss and DiVincenzo 
1998), or with an out-of-plane magnetic field or weak in-plane electric field 
(Burkard et al. 1999). More recently, experiments have controlled J by varying 
the back-gate voltages on two neighboring quantum dots through a large pa- 
rameter regime, independently (Petta et al. 2005a). Here we discuss this last 
method to control J, which has been analyzed in several recent papers (Petta 
et al. 20056, Coish and Loss 2005, Taylor et al. 2006, Stopa and Marcus 2006). 

We consider a double quantum dot in the region of the charge stability 
diagram indicated in the lower inset of Figure Specifically, we consider the 
regime of gate voltages where the double dot contains N = 2 electrons near the 
degeneracy point of the (1,1) and (0,2) charge states, and aim to diagonalizc 
the Hamiltonian Hq + Ht in the basis of three spin triplets and two relevant 
singlets: 



|5(0,2)) 


= 4;4t l vac -> ' 


15(1,1)) 


- _L (j $ _ 

_ V2 V 2i 1T 


lib) 


= \71 (4j.4t + 


\T+) 


= dj-f^it l vac -) i 


|T_) 


= 4i4i l vac -) ■ 



|vac.) , 



(14) 
(15) 

(16) 

(17) 
(18) 



In the absence of additional spin-dependent terms, the triplets are degenerate, 
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with energy i?Tripict = -^(1,1) = — v2A', whereas the two singlet states have 
energies and associated eigenvectors 

kinglet = £TViplet-^(e'±yV) 2 +4i? 2 ), (19) 

^ngict) = cosf$)|S(l,l))+sinf^|5(0,2)), (20) 



-(f) ■ 

Here, A' and e 1 are related to the previous coordinates (A, e) through a simple 
translation of the origin: 



t' ) + U'-U 



(22) 



This gives rise to the Heisenberg exchange for large negative e' (from Equation 
©): 

J(e') = triplet - E+ nglct « ^-^ e'<0, |e'|»2i 12 . (23) 

By pulsing e' = e'(t), the exchange J(e'(t)) can be pulsed on and off again in 
order to implement the vSWAP operation, as described in Section El (see the 
inset of Figure EJ- This operation has now been achieved experimentally with a 
gating time on the order of 180 ps (Petta et al. 2005a), in good agreement with 
the predictions in (Burkard et al. 1999) for an achievable switching time. 



4.4 Initialization of two-spin encoded qubits 

Fluctuations in a nuclear spin environment can lead to rapid decoherence of 
single-electron spin states due to the contact hyperfine interaction (see Sec- 
tion [5J below). The effects of these fluctuations can be reduced, in part, by 
considering a qubit encoded in two-electron singlet |0) = \S(1, 1)) and triplet 
states |1) = \Tq), as defined in equations i|15|) . (|16fl . With this encoding scheme, 
the qubit energy splitting would be provided through the exchange coupling 
(Equation HSJ), and single-qubit rotations could be performed using an inho- 
mogeneous magnetic field (Levy 2002). Two-qubit operations in this scheme 
could be performed, for example, using capacitive coupling due to the relative 
charge distributions of the triplet and singlet states in neighboring double-dots 
(Taylor et al. 2005), although the difference in these charge distributions can 
lead to additional dephasing due to fluctuations in the electrical environment 
(Coish and Loss 2005) (see also Sectional below). An alternative scheme to 
couple such encoded qubits over long distances with optical cavity modes has 
also been proposed (Burkard and Imamoglu 2006). One additional advantage of 
the two-spin encoded qubit scheme is that adiabatic tuning of the gate voltages 
can be used to initialize and readout information stored in the singlet-triplet 
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Figure 7: Leakage (the occupation probability of the state \S(0, 2)) at the end of 
the sweep) due to nonadiabatic transitions after sweeping from e' = 0.7 me V to 
e' = —6 meV. Leakage is given as a function of the characeristic sweep time r sw , 
where e'(t) = eo — ^ tanh (2t/T sw ), with e = — 2.65meV and Ae = 6.7meV. 
We show results for t\i = 5, 7, 10, 20 u.eV. Solid lines show the results of numer- 
ical integration of the Schrodinger equation and dashed lines give exponential 
fits, which decay with the time constants ret, given in Tabled 



basis (Johnson et al. 2005 a, Petta et al. 2005 a). We discuss this initialization 
scheme in the rest of this section. 

We consider the singlet ground-state -Egi nglot ^), given by Equations lf2*UI) and 
l|21|). For large positive detuning, |e'| ^> ti 2 , e' > 0, the mixing angle in Equation 
<|21l) is df_ ss 7r and the singlet ground state is approximately given by \S(Q, 2)). 
For large negative detuning |e'| ^> ii2,e' < 0, we find 6+ sw and the lowest- 
energy singlet is instead given by 15(1, 1)) (see Figure EJ). If the two-electron 
system is allowed to relax to its ground state \S(0, 2)) at large positive detuning 
e and the detuning is then varied adiabatically slowly to large negative values, 
the encoded qubit can be initialized to the state |0) = \S(1, 1)) (see Figure O 
and insets). It is a straightforward exercise to estimate the error in such an 
operation for a two-dimensional Hamiltonian. 

For a linear ramp of e' over an infinite interval (i.e., e' = (Ae/r sw )t, t = 
— oo . . . oo with characteristic switching time r sw to sweep over an interval Ae), 
the result of Zener, for the non-adiabatic Landau-Zener transition probability 
is(Zener 1932) 

f r sw \ hV2Ae 
^ = cxp(-— ), t lz = . (24) 

c 12 



TLZ J 1 47rf? 



Here, the Landau-Zener tunneling probability is controlled in terms of two time 
scales: the switching time r sw for a typical range of Ae and the Landau-Zener 
time constant tlz, which has a strong dependence (oc l/t\^) on the inter-dot 
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h2 (fieV) 


tlz (ns) 


Tfit (ns) 


5 


20. 


7.4 


7 


10. 


3.8 


10 


4.9 


1.9 


20 


1.2 


0.43 



Table 1: Landau-Zener time constant tlz for a linear ramp of e' and the time 
constant Tfit for fits to numerically evaluated data at various values of the tunnel 
coupling ti2- 

tunnel coupling. For a realistic voltage pulse, e' is swept over a hnitc inter- 
val, and the pulse shape, in general, will not be linear for the entire sweep. 
Performing an analysis similar to that used for single-spin gates (Schliemann 
et al. 2001, Requist et al. 2005) for this case, one can perform a numerical inte- 
gration of the time-dependent Schrodinger equation in the subspace formed by 
the two singlets for an arbitrary pulse shape. We have done this for a pulse of the 
form e'(t) = e — tanh (2i/r sw ) , t = — 5t sw . . . 5^^, where we find an approx- 
imately exponential dependence of P on the switching time r sw (see Figure 0). 
Fitting to this exponential dependence, we find a time constant rgt analogous 
to the Landau-Zener time tlz- The time constants tlz (from Equation 11241) ) 
and Tfit from the numerical data in Figure[7]are compared in Table [2 for various 
values of the inter-dot tunnel coupling t%2. The results of Figure Hand Tabled 
suggest that (for this set of parameters) adiabatic switching for initialization or 
readout on a time scale of t sw ^ 1 ns can only be performed without significant 
error if the tunnel coupling t±2 is made larger than t±2 > 20 fieV. It is impor- 
tant to note that this analysis ignores additional effects due to magnetic field 
inhomogeneities, spin-orbit coupling, or the hyperfine interaction, all of which 
can lead to additional singlet-triplet anticrossings (see Section I5"2l below) and 
hence, to additional initialization or readout errors. 

5 DECOHERENCE 

Decoherence is the process by which information stored in a quantum bit is lost. 
There are two time scales used to describe decoherence processes for a spin that 
decays exponentially in the presence of an applied magnetic field. T\ is the 
longitudinal spin decay time, or spin-flip time, which describes the time scale 
for random spin flips: \\) — > | j). 7b , the transverse spin decay time, describes 
the decay of a superposition state a||) + b\l). Both of these time scales are 
important for quantum computing, since both effects lead to qubit errors. 

An experiment performed on an ensemble of systems with different environ- 

1 Note that this type of pulse will generally lead to a smaller value of P for a given set of 
parameters since here, de' /dt < Ae/T sw , whereas for the linear pulse de' /dt = Ae/r Bvi for the 
entire sweep. 
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ments can lead to additional decoherence, beyond that described by the "intrin- 
sic" T 2 time (Slichter 1980). For such an experiment, the ensemble-averaged 
transverse spin decay time is therefore often denoted T 2 * to distinguish it from 
the single-spin decay time. Other symbols such as r c (the correlation time) and 
Tm (the magnetization envelope decay time) are often used to distinguish decay 
that is non-exponential. 

For a quantum-dot-confined electron spin state to decay, it is necessary for 
the spin to couple in some way to fluctuations in the environment. There are 
two important sources of this coupling for electron spins in quantum dots. First, 
the spin-orbit interaction couples electron spin states to their orbital states, and 
therefore makes spins indirectly sensitive to fluctuations in the electric environ- 
ment. Second, the Fermi contact hypcrfine interaction between electrons and 
surrounding nuclear spins in the host material can lead to rapid decay if fluctu- 
ations in the nuclear spin environment are not properly controlled. In the rest 
of this section we discuss recent progress in understanding decoherence due to 
these two coupling mechanisms. 

5.1 Spin-orbit interaction 

For a 2DEG formed in GaAs, the spin-orbit interaction is given in terms of two 
terms: 

H so = a (p x a v - p y (7 x ) + (3 (pyCTy - p x a x ) + O (l P 1 3 ) , (25) 

where a x , y are Pauli matrices and p = {p x ,p y ) is the momentum operator 
in the plane of the 2DEG. The first term, proportional to a, is the Rashba 
(or structure-inversion-asymmetry) spin-orbit coupling term. The Rashba term 
is due to asymmetry in the confining potential and can therefore be tuned 
to some degree with applied gates. The second term, proportional to 0, is 
the Dresselhaus (bulk-inversion-asymmetry) term, and is due to the fact that 
GaAs, which has a zincblende lattice, has no center of inversion symmetry. 
Corrections to this spin-orbit Hamiltonian of order |p| 3 are smaller than the 
linear-momentum terms in quantum dots by the ratio of z-confinement length 
to the quantum-dot Bohr radius, and are negligible in the two-dimensional limit 
(Cerletti et al. 2005). 

Hgo obeys time-reversal symmetry. Thus, in the absence of a magnetic field, 
the ground state of a single electron confined to a quantum dot is twofold de- 
generate due to Kramer's theorem, and H$o alone can not cause decoherence. 
The character of the ground-state doublet does change, however, due to the 
presence of Hso, mixing orbital and spin states. Thus, any fluctuations that 
couple to the orbital degree of freedom can cause decoherence in combination 
with spin-orbit coupling. These fluctuations can come from lattice phonons, 
surrounding gates, electron-hole pair excitations, etc. (Golovach et al. 2004). 
The longitudinal-spin relaxation rate 1/Ti due to spin-orbit coupling and lattice 
phonons has been calculated, and shows a strong suppression for confined elec- 
trons (with large level spacing tkv ) in weak magnetic fields: 1/Ti cx B 5 / (hu! ) 4 
(Khaetskii and Nazarov 2000, Khaetskii and Nazarov 2001). This calculation 
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has been extended to a larger range of magnetic fields, showing that the B- 
ficld dependence of 1/Ti saturates and is then suppressed when the phonon 
wavelength is comparable to the dot size (Golovach et al. 2004). Further, this 
calculation has also been extended to include the transverse spin decay time 
due to spin-orbit interaction alone, showing that dephasing is limited by relax- 
ation, or Ti — 2T\ to leading order in the spin-orbit coupling, independent of 
the particular source of fluctuations (Golovach et al. 2004). Additionally, 1/Ti 
has been shown to have a strong dependence on the magnetic field direction, 
relative to the crystal axes (Fal'ko et al. 2005), shows a strong enhancement 
near avoided level crossings, which may allow independent measurements of the 
Rashba and Dresselhaus coupling constants (Bulaev and Loss 2005 a), and plays 
a role in phonon-assisted cotunneling current through quantum dots (Lehmann 
and Loss 2006). The relaxation rate of quantum-dot-confined hole spins due 
to spin-orbit coupling and phonons has also been investigated. In some cases, 
recent work has shown that the hole spin relaxation time may even exceed the 
relaxation time of electron spins (Bulaev and Loss 2005 b). In addition to lattice 
phonons, electric field fluctuations can result from the noise in a QPC readout 
device, which results in spin decoherence when considered in combination with 
spin-orbit coupling. This mechanism shows a strong dependence of the decoher- 
ence rate ~ 1/r 6 on the dot-QPC separation r (Borhani et al. 2005), and can 
therefore be controlled with careful positioning of the readout device. 

Measurements of relaxation times for single electron spins have been per- 
formed in gated lateral quantum dots (Hanson et al. 2003, Elzerman et al. 2004), 
giving a T\ time in good agreement with the theory of ref. (Golovach et al. 2004) 
and in self-assembled quantum dots (Kroutvar et al. 2004), which confirmed the 
expected magnetic field dependence: 1/Ti oc B 5 (Khaetskii and Nazarov 2001). 
Additionally, singlet-triplet decay has been measured in single vertical (Fujisawa 
et al. 2002), and lateral (Hanson et al. 2005) dots, as well as lateral double dots 
(Petta et al. 20056, Johnson et al. 2005a). 

There is a general consensus that spin relaxation for quantum-dot-confined 
electrons proceeds through the spin-orbit interaction and phonon emission at 
high magnetic fields. However, in weak magnetic fields, and for the transverse 
spin decay time T 2 , there are stronger effects in GaAs. These effects are due to 
the contact hyperfine interaction between confined electron spins and nuclear 
spins in the surrounding lattice. 

5.2 Hyperfine interaction 

For a collection of electrons in the presence of nuclear spins, the Fermi contact 
hyperfine interaction reads 

H hf = AvJ2U-S(r k ), S(r fc ) = i ^ V>£(rfc)<wVv(r fe ), (26) 

fc <7,<t' = 1,1 

where A is the hyperfine coupling strength, v is the volume of a crystal unit cell 
containing one nuclear spin, 1^ is the spin operator for the nuclear spin at site 
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k, S(rfe) is the electron spin density at the nuclear site, given in terms of field 

operators ipo-i?) that satisfy the anticommutation relations |^v(r), ip^, (r')| = 

S(r — T')S^ a i, {?/v(r), f/v (r')} = 0, and we have denoted matrix elements by 
o a o' = (cr| cr |cr') , where a = {a x ,a y ,a z ) is the vector of Pauli matrices. The 
significance of the general form given in Equation Ij26(l is that there can be 
an interplay of orbital and spin degrees of freedom due to the contact hyperfine 
interaction. When the orbital level spacing is not too large, this interplay can be 
the limiting cause of electron-spin relaxation (Erlingsson et al. 2001, Erlingsson 
and Nazarov 2002) in weak magnetic fields, where the spin-orbit interaction is 
less effective, and leads to enhanced nuclear spin relaxation in the vicinity of 
sequential-tunneling peaks for a quantum dot connected to leads, where S(r) 
fluctuates significantly (Lyanda-Geller et al. 2002, Huttel et al. 2004). 

The orbital level spacing in lateral quantum dots is usually much larger than 
the typical energy scale of Hm . In this case, it is possible to solve for the orbital 
envelope wavefunction \I/o(r) i n the absence of the hyperfine interaction, and 
write an effective hyperfine Hamiltonian for a single electron confined to the 
quantum-dot orbital ground state: 

ff M ,o«b.o-So, h = ^^|*o(r fc )| 2 I fc , (27) 

where here So is the spin-1/2 operator for a single electron in the quantum-dot 
orbital ground state. The primary material used to make lateral quantum dots 
is GaAs. All natural isotopes of Ga and As carry nuclear spin / = 3/2. Each 
isotope has a distinct hyperfine coupling constant, but the average coupling 
constant, weighted by the relative abundance of each isotope in GaAs gives 
A w 90 fieV (Paget et al. 1977). 

Dynamics under Hm .o have now been studied extensively under many various 
approximations and in many parameter regimes. Here we give a brief account of 
some part of this study. For an extensive overview, see reviews in (Schliemann 
et al. 2003, Cerletti et al. 2005). The first analysis of the influence of Equation 
l|27|) on quantum-dot electron spin dynamics showed that the long-time longi- 
tudinal spin- flip probability, P n w l/p 2 N (Burkard et al. 1999) was suppressed 
in the limit of large nuclear spin polarization p and number of nuclear spins in 
the dot, N. Subsequently, an exact solution for the case of a fully-polarized 
nuclear spin system (p = 1) has shown that both the longitudinal and trans- 
verse components of the electron spin decay by a fraction ~ 1/N according 
to a long-time power law ~ 1/t 3 / 2 on a time scale of r ~ hN/A (Khaetskii 
et al. 2002) (hN/A ~ 1/is for a GaAs dot containing N ~ 10 5 nuclei). This 
exact solution for p = 1, which shows a non-exponential decay, demonstrates 
that the electron spin decay is manifestly non-Markovian since the time scale 
for motion in the nuclear-spin bath is much longer than the decay time scale of 
the electron spin. For unpolarized systems, the ensemble averaged mean-field 
dynamics show a transverse spin decay on a time scale r ~ h^/N /A ~ 5 ns 
(Khaetskii et al. 2002, Merkulov et al. 2002). The exact solution has been ex- 
tended to the case of nonzero polarization p ^ 1 using a generalized master 
equation, valid in the limit of large magnetic field or polarization p 3> 1/VA" 
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(Coish and Loss 2004). This work has shown that, while the longitudinal spin 
decay is bounded by ~ l/p 2 N, due to the quantum nature of the nuclear field, 
the transverse components of spin will decay to zero in a time t c fts 5 ns/yl — p 2 
(without ensemble averaging and without making a mean-field ansatz) , unless an 
electron spin echo sequence is performed or the nuclei are prepared in an eigen- 
state of the operator h$ through measurement (Coish and Loss 2004). There are 
several recent suggestions for methods that could be used to measure the oper- 
ator h,Q (Giedke et al. 2005, Klauser et al. 2005, Stepanenko et al. 2005) in order 
to extend electron spin decoherence. Once the nuclear spin system is forced 
into an eigenstate of Hq, the lowest-order corrections for large magnetic field 
still show incomplete decay for the transverse spin (Coish and Loss 2004), sug- 
gesting that dynamics induced by the nuclear dipolar interaction may limit spin 
coherence in this regime (de Sousa and Das Sarma 2003), although higher-order 
corrections have been reported to lead to complete decay (Deng and Hu 2005), 
even when the nuclear spin system is static. There have been several ef- 
forts to understand the hyperfine decoherence problem numerically (Schliemann 
et al. 2002, Shenvi et al. 2005&), and other studies have investigated electron 
spin-echo envelope decoherence under the hyperfine interaction alone (Coish and 
Loss 2004, Shenvi et al. 20056, Shenvi et al. 2005a) or the combined influence 
of hyperfine and nuclear dipolar interactions (de Sousa and Das Sarma 2003, de 
Sousa et al. 2005, Witzel et al. 2005, Yao et al. 2005, Yao et al. 2006). Other 
approaches to understanding the hyperfine decoherence problem include semi- 
classical theories that replace the quantum nuclear field by a classical dynamical 
vector (Erlingsson and Nazarov 2004, Yuzbashyan et al. 2004) or a classical dis- 
tribution function (Al-Hassanieh et al. 2005). 

Experiments on electron spin decoherence in single quantum dots (Bracker 
et al. 2005, Dutt et al. 2005) and double quantum dots (Petta et al. 2005a, 
Koppens et al. 2005) have now confirmed that the ensemble-averaged electron 
spin dephasing time is indeed given by r ~ h\f~N j A ~ 10 ns. 

For two electron spins confined to a double quantum dot, the hyperfine 
Hamiltonian (Equation J2HJ)) can be cast in the form (Coish and Loss 2005) 

ffhf,dd = e,5f +5^h,-S,; S i = i £ d\ a a^d l<y ,, (28) 

l ff,<r'=T4 

= e z S 7 { + S • h + SS ■ Sh, (29) 

where e 2 = g/iBB is the Zeeman splitting, d 1 ( 2 ) cr annihilates an electron in the 

single-particle orbital state with envelope wavefunction ^ 1 ( 2 )(r) and spin a, we 

define h = (hi + h 2 ) /2, 5h = (hi + h 2 ) /2, where the quantum nuclear field 

■ 1 2 

operators are hi(2) = Av^2 k hJji(2)(rfc) the sum of electron spins is S — 
Si + S 2 and the difference is SS = Si — S2. While the sum S conserves the total 
squared electron spin, and can only couple states of different z-projection (e.g. 
|To) to |T±)), the difference <5S does not preserve the total spin, and therefore 
couples singlet to triplet (e.g. |5(1,1)) to |Tq) and \T±)) and will therefore 
lead to anticrossings in the energy level spectrum, where \S(1, 1)) and \T±) 
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or |Xb) cross. Adding Equation (|29J) to the previous double-dot hamiltonian, 
Hdd = He + Ht + -ffhf.dd, and making a mean- field ansatz for the nuclear field 
operators, i.e., replacing operators by their expectation values: h — > (h), 2 leads 
to the energy level spectrum shown in Figure In the limit of large Zeeman 
splitting e z and large negative detuning e', an effective two- level Hamiltonian can 
be derived in the subspace of lowest-energy singlet and S z = triplet (\S) , |Tb)) 
(Coish and Loss 2005): 



An exact solution can be found for pseudospin dynamics in the two-dimensional 
subspace of \S) and \Tq) under the action of -ffdd,cff- This solution shows that a 
singlet-triplet correlator undergoes an interesting power-law decay in a charac- 
teristic time scale that can be extended by increasing J (Coish and Loss 2005), 
and has been verified in experiment (Laird et al. 2005). As is true for the 
transverse components of a single electron spin, the singlet-triplet correlator 
shows a rapid decay if the nuclear spin environment is not in an eigenstate of 
the relevant nuclear field operator (in this case, Sh z ). The decay time can be 
significantly extended by narrowing the distribution in 5h z eigenstates through 
measurement (Klauser et al. 2005) or by performing a spin-echo sequence (Petta 
et al. 2005a). Remaining sources of dephasing include the corrections to -ffdd.cff 
(of order 1 / ' e z , which can not be removed easily) and fluctuations in the electro- 
static environment, although the effect of these fluctuations can be removed to 
leading order at zero-derivative points for the exchange interaction (Coish and 
Loss 2005), where: 



Recent calculations suggest that these zero-derivative points should be achiev- 
able with appropriate control of the confinement potential or magnetic field (Hu 
and Das Sarma 2006, Stopa and Marcus 2006). 

Since the hyperfine interaction does not preserve the total spin quantum 
number of electrons, this interaction plays a very important role in studies on 
spin-dependent transport. In particular, spin blockade (Weinmann et al. 1995, 
Weinmann 2003) occurs in double quantum dots (Ono et al. 2002) when tunnel- 
ing is allowed betwen spin-singlets 1)) — > \S(0, 2)), but not between spin 
triplets |T(1, 1)) -y* \T(0, 2)), because of a large energy cost due to orbital level 
spacing and the Pauli principle. This blockade allows for the extraction of fea- 
tures at energy scales much less than temperature, making it an ideal param- 
eter regime in which to perform spectroscopy on double dots (Pioro-Ladriere 
et al. 2003) and spin-resonance experiments, which previously suffered from 
"heating" effects in single dots (van der Wiel et al. 2003, Koppens et al. 2006). 
The hyperfine interaction mixes the 1)) and |T(1, 1)) states, allowing trans- 
port, and effectively removing spin blockade when these states are nearly de- 
generate. This behavior leads to a number of intriguing effects, including stable 

2 In general, great care should be taken in making such a replacement. See the discussion, 
for example, in (Coish and Loss 2005). 
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undrivcn oscillations in transport current (Ono and Tarucha 2004, Erlingsson 
et al. 2005), and a striking magnetic-field dependence of leakage current, which 
allows the extraction of information about the nuclear spin system (Koppens 
et al. 2005, Jouravlev and Nazarov 2006). Even-odd effects in the spin blockade 
of many-electron quantum dots have further revealed the shell-filling illustrated 
in Figure 01 (Johnson et al. 20056). 

The influence of spin-dependent terms, causing decohcrcnce or unwanted 
evolution, is a central issue in quantum-dot spin quantum computing. The re- 
quirements for fault-tolerant quantum information processing are very stringent. 
This raises the bar for required understanding of these environmental influences 
to a very high level, and guarantees that quantum-dot spin decoherence will 
remain a challenge for some time to come. 

6 ENTANGLEMENT GENERATION, DISTIL- 
LATION, AND DETECTION 

In addition to the usual requirements for control and coherence, to demon- 
strate the true quantum nature of qubits, there have been many suggestions 
to create and measure nonlocal multiparticle entanglement of electron spins 
in nanostructures (DiVincenzo and Loss 1999, Burkard et al. 2000, Loss and 
Sukhorukov 2000, Choi et al. 2000, Egues et al. 2002, Burkard and Loss 2003, 
Samuelsson et al. 2004, Recher et al. 2001, Lesovik et al. 2001, Melin 2001, 
Costa and Bose 2001, Oliver et al. 2002, Bose and Home 2002, Recher and 
Loss 2002, Bena et al. 2002, Saraga and Loss 2003, Bouchiat et al. 2003, Recher 
and Loss 2003, Beenakker and Schoenenberger 2003, Saraga et al. 2004, Egues 
et al. 2005). These proposals include suggestions to extract spin singlets from 
a superconductor through two quantum dots (Recher et al. 2001) or nanotubes 
(Recher and Loss 2002, Bena et al. 2002), or to create entanglement near a mag- 
netic impurity (Costa and Bose 2001), through a single quantum dot (Oliver 
et al. 2002), from biexcitons in double quantum dots (Gywat et al. 2002), or 
through a triple dot (Saraga and Loss 2003). It may also be possible to dis- 
till entanglement (Bennett et al. 1997) from an unentangled Fermi gas through 
Coulomb scattering in a 2DEG (Saraga et al. 2004). 

As well as providing a proof of quantum mechanical behavior, entanglement 
can be used as a resource for measurement-based quantum computing. Some 
measurement-based schemes rely on the creation of highly-entangled cluster 
states (Raussendorf and Briegel 2001), which could be generated in quantum- 
dot arrays using the Heisenberg exchange interaction (Borhani and Loss 2005). 
Other measurment-based schemes generate entanglement through partial Bell- 
state (parity) measurements (Beenakker et al. 2004), which could also be imple- 
mented for spins in quantum dots using spin-to-charge conversion (Engcl and 
Loss 2005). Independent of the method used, the generation or purification 
and subsequent detection of entangled electron spins would present a significant 
milestone on the road to a working quantum-dot quantum computer. 
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7 CONCLUSIONS AND OUTLOOK 

We have presented some of the theoretical and experimental challenges to quantum- 
dot quantum computing with electron spins. The last few years have seen an 
extremely rapid rate of progress in experiments which show that many of the 
required elements of a spin-based quantum-dot quantum computer can be re- 
alized in principle. The most significant advances include the reduction of the 
number of electrons confined to gated quantum dots down to a single electron 
(Ciorga et al. 2000), the demonstration (Hanson et al. 2003) and improvement 
(Hanson et al. 2005) of electron spin readout in gated lateral dots, which has led 
to the measurement of a spin T\ time (Elzerman et al. 2004), the demonstration 
of the Vswap operation, allowing for the extraction of an ensemble-averaged 
T| time, and spin-echo methods to extend the decay time within a two-spin 
encoded subspace (Petta et al. 2005 a), and most recently the demonstration of 
single-spin rotations under resonant conditions (Koppens et al. 2006). 

To demonstrate viability of the Loss-DiVincenzo proposal, more experiments 
are needed. Although the Loss-DiVincenzo proposal is scalable in principle, it 
remains to be seen if there are significant practical obstacles to scaling-up the 
number of electrons involved well beyond two. 
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